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A  laboratory  experiment  was  conducted  to  measure  the  speed  of  sound  in  an  artificial 
water-saturated  granular  sediment  composed  of  cleaned  and  sorted  medium-grained  sand  and 
degassed  distilled  water.  The  experiment  was  conducted  within  a  range  of  frequencies  where 
dispersion  is  predicted  by  a  number  of  existing  models.  Between  2  and  4  kHz,  the  sound  speed  was 
inferred  from  measurements  of  the  resonance  frequencies  of  a  thin-walled  cylindrical  container 
filled  with  the  material.  An  elastic  waveguide  model  was  used  to  account  for  the  effect  of  the  finite 
impedance  of  the  walls,  although  this  effect  was  found  to  be  small.  From  20  to  300  kHz,  the  sound 
speed  was  obtained  directly  from  time-of-flight  measurements  within  the  sediment.  Dispersion  in 
close  agreement  with  the  Williams  effective  density  fluid  model  [K.  L.  Williams,  J.  Acoust.  Soc. 
Am.  110,  2276-2281  (2001)]  was  observed.  €)  2007  Acoustical  Society  of  America. 
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I.  INTRODUCTION 

A  complete  understanding  of  sound  propagation  in 
water-saturated  granular  sediments  has  not  been  achieved.  A 
number  of  experimental  studies  have  appeared  in  the  litera¬ 
ture  that  report  either  direct  or  indirect  evidence  of  sound 
speed  dispersion  in  this  medium.1  6  Such  results  generally 
support  the  Biot-Stoll  model  (and  its  variants)7-12  for  sound 
propagation.  The  Biot-derived  models  predict  a  low- 
frequency  sound  speed,  followed  by  an  increase  to  a  higher 
sound  speed  beyond  some  transition  frequency.  There  are 
also  a  number  of  experimental  studies17 14  that  found  no  evi¬ 
dence  of  sound  speed  dispersion  and  support  Hamilton’s  dis¬ 
persionless  fluid  model.1  The  dispersion  predicted  by  the 
Biot  models  occurs  over  a  range  of  about  two  orders  of  mag¬ 
nitude  in  frequency  (from  around  100  Hz  to  10  kHz  for 
medium-grained  sand,  for  example).  It  is  experimentally  dif¬ 
ficult  to  obtain  measurements  on  a  uniform  body  of  sediment 
using  a  single  technique  that  spans  the  entire  range  of  inter¬ 
est.  Further,  the  predicted  dispersion  is  on  the  order  of  10% 
and  achieving  experimental  accuracy  significantly  greater 
than  this  has  been  difficult.  Surface  and  volume  inhomoge¬ 
neities  in  natural  sediments  also  cause  uncertainty  in  experi¬ 
mental  measurements.4,16 

A  series  of  laboratory  experiments  designed  to  overcome 
some  of  these  difficulties  were  performed  on  a  water- 
saturated  sand  sediment.  The  sediment  was  contained  within 
a  thin-walled  cylindrical  vessel  large  enough  to  permit  time- 
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of-flight  sound  speed  measurements  within  the  bulk  sediment 
above  20  kHz.  Lower  frequency  sound  speeds  were  inferred 
from  the  measured  resonance  frequencies  of  the  cylindrical 
sediment-filled  container.  For  this  analysis,  transverse  waves 
within  the  sediment  were  ignored.  The  sediment  was  mod¬ 
eled  as  a  fluid,  but  the  sediment  sound  speed  was  not  con¬ 
strained  to  be  constant.  The  walls  of  the  cylinder  were  ini¬ 
tially  approximated  as  a  pressure  release  boundary  (details  of 
the  assumption’s  validity  are  discussed  in  the  following).  The 
cylinder  was  placed  upon  a  layer  of  insulation  foam,  which 
yielded  a  pressure  release  lower  boundary  condition.  The 
water/air  interface  provided  a  pressure  release  boundary  at 
the  top  of  the  sediment.  Therefore,  a  pressure  release  wave¬ 
guide  model  was  used  to  relate  the  measured  resonance  fre¬ 
quencies  to  the  effective  sediment  sound  speed.  In  order  to 
identify  individual  modes,  and  hence  infer  frequency- 
dependent  sound  speeds  from  resonances,  both  symmetric 
and  asymmetric  excitation  of  the  sediment  cylinder  was 
used.  The  spatially  dependent  phase  of  the  acoustic  pressure, 
relative  to  the  drive  signal,  was  employed  to  aid  in  mode 
identification.  Dispersion  was  observed  and  these  results  are 
presented  and  compared  to  the  effective  density  fluid 
model,1  hereafter  referred  to  as  EDFM.  The  effect  of  the 
finite  impedance  imparted  by  the  cylindrical  wall  of  the  sedi¬ 
ment  container  is  also  discussed.  A  model  of  an  elastic- 
walled  waveguide  was  used  to  relate  the  sediment  sound 
speeds  that  were  inferred  from  axisymmetric  modes  in  the 
tank  to  the  values  that  would  be  observed  in  an  unconfined 
free  field.  The  effect  was  found  to  be  small  (less  than  0.53%) 
at  the  lowest  experimental  frequency  and  decreased  as  the 
frequency  increased. 
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FIG.  1 .  A  schematic  diagram  of  the  time-of-flight  measurement  apparatus. 

II.  DESCRIPTION  OF  THE  EXPERIMENT 

A.  Sediment  container 

The  sediment  was  contained  in  a  right-circular- 
cylindncal  tank  possessing  walls  made  of  high-density  poly¬ 
ethylene  (HDPE).  The  mean  inner  diameter  was  0.686  m,  the 
height  of  the  sediment  within  the  vessel  was  0.80  m,  and  the 
mean  wall  thickness  was  0.48  cm.  The  tank,  which  had  a  flat 
bottom,  was  placed  upon  a  10-cm-thick  styrofoam  panel 
which,  in  turn,  rested  on  a  concrete  floor. 

The  sediment  was  prepared  in  the  following  manner.  Fil¬ 
tered  fresh  water  was  added  to  the  tank  and  degassed  in  situ 
by  heating  it  to  60  °C  using  several  immersion  heaters.  Cir¬ 
culation  was  imposed  during  this  process  with  a  motorized 
stirring  impeller.  The  water  was  kept  at  60  °C  for  several 
hours  during  which  time  the  dissolved  gas  concentration  at¬ 
tained  equilibrium.  Then,  approximately  500  kg  of  a  com¬ 
mercial  silica  blasting  sand  was  slowly  added  to  the  water,  as 
the  stirring  continued.  The  sediment  was  then  allowed  to 
cool  to  room  temperature,  over  a  period  of  24  h.  As  the 
system  cooled,  the  water  phase  became  undersaturated  with 
gas  and  thus  drew  into  solution  any  gas  pockets  entrained  by 
the  introduction  of  the  sand,  resulting  in  a  fully  saturated 
sediment.  A  sieve  analysis  conducted  on  approximately 
500  ml  of  the  sand  yielded  grain  sizes  ranging  from  150  to 
600  /xmy  with  a  mean  size  of  300  /im. 

B.  High-frequency  time-of-flight  measurements 

The  sound  speed  within  the  sediment  was  measured  us¬ 
ing  a  time-of-flight  technique  at  20,  50,  100,  and  300  kHz.  A 
schematic  diagram  of  this  experiment  is  shown  in  Fig,  1 .  The 
signal  source  was  a  function  generator  programmed  to  pro¬ 
duce  sinusoidal  pulses  of  various  length,  depending  upon  the 
frequency,  to  avoid  multipath  interference.  The  source  signal 
was  bandpass  filtered  (24  dB/octave)  at  the  pulse  center  fre¬ 
quency  with  a  Krohn-Hite  model  34A  filter  and  then  ampli¬ 
fied  by  a  Krohn-Hite  model  7500  wide-band  power  amplifier 
and  directed  to  a  Briiel  &  Kjaer  model  8103  miniature  hydro¬ 
phone,  operating  as  a  source  (labeled  S  in  Fig.  1).  Two  more 
8103  hydrophones  were  positioned  10.0  and  17.6  cm  away 
from  the  source  on  a  horizontal  line  (parallel  to  the  sediment 
surface),  at  a  depth  of  15.0  cm  The  latter  two  hydrophones 
acted  as  receivers  (labeled  R1  and  R2  in  Fig.  1).  The  relative 
positions  of  the  source  and  receivers  were  constrained  to 
within  ±0.25  mm  by  mounting  each  hydrophone  and  its 
cable  inside  a  stiff  stainless  steel  tube.  These  tubes  were 


FIG.  2,  A  schematic  diagram  of  the  resonance  frequency  apparatus. 

mounted  to  a  steel  frame  and  fastened  securely  in  place.  The 
R1-R2  separation  distance  was  set  and  measured  with  a  mi¬ 
crometer  caliper  prior  to  inserting  the  hydrophones  into  the 
sediment.  The  rigidity  of  the  mounting  scheme  was  sufficient 
to  prevent  caliper-measureable  deflection  during  the  inserting 
process.  The  hydrophone-enclosing  tubes  were  filled  with 
degassed  water  to  minimize  acoustic  disturbance  of  the  sedi¬ 
ment. 

The  signals  from  R1  and  R2  were  conditioned  with  a 
Brtiel  &  Kjaer  model  2692  charge  amplifier,  then  bandpass 
filtered  (24  dB/octave)  at  the  pulse  center  frequency  with 
Krohn-Hite  model  34A  filters.  The  filtered  source  and  re¬ 
ceive  signals  were  acquired  with  a  digital  oscilloscope.  At 
each  frequency,  100  pings  were  collected  and  the  averaged 
signals  were  sent  to  a  computer  for  processing.  The  R1  and 
R2  signals  were  cross-correlated  to  obtain  the  time-of-flight 
and  then  sound  speed  was  calculated  from  the  known  sensor 
separation  distance.  The  uncertainty  in  sound  speed  was  cal¬ 
culated  from  the  previously  stated  position  uncertainty,  and 
the  digitization  sample  rate  and  found  to  be  between  ±6.6 
and  ±9.9  m/s  at  300  and  20  kHz,  respectively. 

C.  Low-frequency  resonance  measurements 

The  apparatus  for  this  portion  of  the  experiment  is 
shown  in  Fig.  2.  A  custom-fabricated  Tonpilz-type  transducer 
was  used  as  a  source.  It  is  described  in  Ref.  17  and  has  a 
circular-shaped  radiating  face  with  an  accelerometer  imbed¬ 
ded  in  it,  and  therefore  can  be  driven  at  constant  velocity 
amplitude  by  suitable  monitoring  of  the  accelerometer  signal 
and  modification  of  the  excitation  voltage  amplitude.  The 
source  was  placed  face-down  on  the  surface  of  the  sediment 
in  one  of  two  positions,  to  be  described  later.  The  source 
signal  was  a  band-limited  (2-5  kHz)  periodic  chirp,  pro¬ 
vided  by  a  Hewlett-Packard  894 10A  vector  signal  analyzer 
(VSA)  and  amplified  by  a  Crown  CE  1000  power  amplifier. 
The  chirp  amplitude  was  modified  onboard  the  VSA  to  pro¬ 
duce  a  source  velocity  amplitude  that  was  uniform  across  the 
frequency  range,  as  discussed  earlier.  The  acoustic  pressure 
was  measured  in  the  sediment  with  B&K  8103  hydrophones 
at  two  positions,  labeled  R1  and  R2  in  Fig.  2.  The  depth  was 
6.7  cm  and  each  sensor  was  positioned  10  cm  from  the  ves¬ 
sel  wall.  When  observed  from  above,  R1-S-R2  formed  a 
straight  line.  The  signals  from  R1  and  R2  were  conditioned 
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TABLE  I.  Specific  acoustic  impedance  of  sediment  and  tank  materials. 


Material 

Specific  acoustic  impedance  (Pas/m) 

Sediment 

3.5  X106 

Tank  wall  material,  HDPE 

2.5  X  106 

Air 

415 

by  the  B&K  2692  charge  amp,  which  included  bandpass  fil¬ 
tering  from  10  to  10  kHz,  and  digitized  by  the  VSA.  The 
complex  spectrum  of  each  hydrophone  signal  (with  phase 
referenced  to  the  excitation  signal)  was  calculated  onboard 
the  analyzer  using  100  spectral  averages  and  converted  to 
complex  acoustic  pressure  spectra  using  the  calibrated  hy¬ 
drophone  sensitivities. 

III.  ACOUSTIC  MODEL  OF  SEDIMENT  AND 
CONTAINER 

An  acoustic  model  of  the  sediment/container  system  is 
needed  to  infer  sound  speeds  from  measured  resonance  fre¬ 
quencies.  In  this  analysis,  the  sediment  was  considered  to  be 
an  effective  fluid  with  wave  number  /c,  which  is  a  valid  ap¬ 
proach  for  granular  sediments.1  We  assume  the  phase  speed 
can  be  frequency  dependent,  but  Im[&]<^Re[£:],  so  that  the 
resonance  frequencies  depend  primarily  upon  the  phase 
speed  and  have  weak  dependence  upon  the  attenuation. 
Transverse  waves  in  the  sediment  were  ignored. 

Initially,  we  ignore  the  container  walls  and  consider  the 
sides  of  the  cylinder  to  have  a  pressure-release  boundary 
condition  under  the  following  justification.  (1)  Typical  values 
of  the  specific  acoustic  impedance  of  the  materials  involved 
are  given  in  Table  I.  Based  on  these  values  and  the  wall 
thickness  (0.48  cm),  the  three-medium  reflection  coefficient 
R ,  given  by  Eq.  6.3.7  in  Ref.  18  for  pressure  waves  incident 
upon  the  wall,  is  /?=  1.00 176.5°  ±  1.5°  between  2  and 
5  kHz.  (2)  This  was  qualitatively  verified  by  the  observation 
of  inverted  reflections  from  the  walls  during  the  time-of- 
flight  measurements,  and  (3)  by  the  absence  of  plane-wave 
longitudinal  mode  resonances  in  the  system.  The  air- 
sediment  interface  at  the  top  of  the  container  is  clearly  pres¬ 
sure  release.  Finally,  the  container  was  placed  upon  a  10-cm- 
thick  styrofoam  building  insulation  panel  to  approximate  a 
pressure  release  boundary  condition  on  the  bottom.  For  the 
cylindrical  coordinate  system  shown  in  Fig.  1  and  sinusoidal 
excitation,  the  acoustic  pressure  in  a  right-circular  cylinder 
with  p- 0  boundary  conditions  at  all  surfaces  is 

Pmnt^r,  6,Z,t)  =  AmnNJm(amnr/a ) 

Xcos(m0- ^0m)sin(Ar7rz//i)^/aj/,  (1) 

and  has  resonance  frequencies 


where  amn  is  the  nth  root  of  the  Bessel  function  Jm(x)y  c  is 
the  effective  speed  of  sound  in  the  sediment,  a  is  the  tank 
radius,  h  is  the  sediment  height,  AmnN  is  the  modal  amplitude 
coefficient,  and  00m  is  a  phase  angle  which  depends  on  the 
source  condition.  Equations  (1)  and  (2)  were  adapted  from 


FIG.  3.  The  resonance  frequencies  for  the  six  lowest-order  mode  families 
for  a  constant  sound  speed  of  1730  m/s.  The  families  are  labeled  m,n 
around  the  top  and  right-hand  side  of  the  plot.  The  plane  wave  mode  would 
be  0,0,7/  if  it  existed. 

Ref.  19.  A  plot  of  resonance  frequencies  versus  longitudinal 
mode  number  N  is  shown  in  Fig.  3  for  the  six  lowest-order 
mode  families.  Five  modes,  in  order  of  ascending  frequency, 
are  labeled  with  roman  numerals  from  I  to  V.  Note  that 
above  mode  III,  multiple  modes  appear  with  similar  eigen- 
frequencies. 

Calculation  of  sound  speed  from  a  measured  spectral 
peak  requires  one  to  identify  the  specific  eigenmode  and  its 
m,n,N  values.  Mode  family  0,1, /V  is  axisymmetric  and 
mode  family  1,1, /V  is  asymmetric  about  the  center  of  the 
sediment  cylinder,  as  shown  in  Fig.  4.  These  two  modes  were 
selectively  excited  by  using  the  source  positions  shown  in 
Fig.  4.  To  excite  mode  family  0, 1  ,N,  the  source  was  placed 
at  the  center  of  the  sediment  cylinder.  To  excite  mode  family 
1,1,  /V,  the  source  was  placed  halfway  between  the  center 
and  the  container  wall.  The  relative  phase  of  the  acoustic 
pressure  observed  at  R1  and  R2  was  used  to  help  identify  the 
mode  associated  with  an  observed  resonance  peak.  For  ex¬ 
ample,  Eq.  (1)  yields  acoustic  pressures  at  R1  and  R2  that  are 
in  phase  for  modes  I,  II,  and  V,  and  180°  out  of  phase  for 
modes  III  and  IV. 

mode  family 

0,1  A  1,1, 


normalized  acoustic  pressure 


FIG.  4.  The  mode  shapes  given  by  Eq.  (1)  for  two  mode  families.  The 
source  positions  for  each  case  are  labeled  S.  The  receiver  positions,  which 
remain  the  same  in  both  cases,  are  labeled  R1  and  R2. 
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resonance 


time-of-flight 


FIG.  5.  Acoustic  pressure  spectra  obtained  at  position  R1  for  both  symmet¬ 
ric  and  asymmetric  excitation.  Spectra  measured  at  position  R2  were  similar. 
The  acoustic  pressure  measured  between  R1  and  R2  was  in  phase  (within 
±1°)  at  the  resonance  peaks  labeled  I,  II,  and  V,  and  was  out  of  phase  (by 
170°)  at  resonances  III  and  IV,  as  indicated  by  ©  and  ©. 


IV.  RESULTS 

A  total  of  four  averaged  spectra  were  obtained  during 
the  resonance  frequency  experiment.  These  consisted  of  the 
acoustic  pressure  measured  at  R1  and  R2,  for  both  of  the 
source  positions  shown  in  Fig.  4.  The  pressure  spectra  ob¬ 
tained  at  R1  are  shown  in  Fig.  5.  For  the  spectral  peaks 
labeled  I,  II,  and  V,  there  was  less  than  one  degree  phase 
difference  between  R1  and  R2.  These  three  spectral  peaks 
are  associated  with  mode  family  0,1  ,/V  and  correspond  with 
the  symmetric  eigenmodes  labeled  I,  II,  and  V  in  Fig.  3. 
Returning  to  Fig.  5,  the  spectral  peaks  labeled  III  and  IV 
exhibited  an  R1-R2  phase  difference  of  170°.  These  peaks 
correspond  to  the  asymmetric  modes  III  and  IV,  shown  in 
Fig.  3  and  associated  with  mode  family  1,1,  A.  Effective 
sediment  sound  speeds  c  were  calculated  from  each  of  the 
identified  peaks  in  Fig.  5  using  Eq.  (2).  Although  not  shown, 
similar  results  were  obtained  from  the  sensor  at  R2. 

Both  the  resonance-based  and  time-of-flight-based 
sounds  speeds  are  shown  in  Fig.  6.  The  speeds  are  normal¬ 
ized  by  the  speed  of  sound  in  distilled  water  at  the  ambient 
temperature  (23.5  ±0.5  °C).  The  sound  speed  predicted  by 
the  effective  density  fluid  model12  is  also  shown  in  Fig.  6. 
The  pertinent  equations  of  the  model  are  shown  in  the  fol¬ 
lowing  for  convenience.  The  sediment  parameters  used  for 
the  model  calculation  are  shown  in  Table  II.  These  were 
primarily  taken  from  Ref.  12  but  the  following  were  obtained 
from  measurements:  The  water  sound  speed  was  calculated 
from  the  experimental  temperature.  The  porosity  of  four 
samples  was  determined  by  the  water  evaporation  method. 
The  permeability  was  measured  using  a  Soiltest  K-605  com¬ 
bination  permeameter  modified  to  accommodate  a  13-cm- 
long,  6.1-cm-diam  core.  Ten  constant-head  permeability 
tests"  were  performed  on  each  of  two  samples. 

V.  SEDIMENT  SOUND  SPEED  FROM  THE  EFFECTIVE 
DENSITY  FLUID  MODEL 

The  effective  density  fluid  model  of  sound  propagation 
in  water-saturated  granular  sediments  is  described  in  detail 
elsewhere,1  but  the  equations  for  prediction  of  the  sediment 


FIG.  6.  Normalized  sound  speed  measurements  are  compared  to  the  EDFM 
(solid  line).  Resonance  measurements  are  indicated  with  circles  and  squares 
(R1  and  R2,  respectively,  with  source  in  center)  and  crosses  and  diamonds 
(R1  and  R2,  respectively,  source  off-center,  near  R2).  Vertical  error  bars 
near  2  and  4  kHz  represent  sound  speed  uncertainty  due  to  uncertainties  in 
the  dimensions  of  the  tank,  but  are  only  shown  for  the  minimum  and  maxi¬ 
mum  sound  speed  measurements,  so  that  the  remaining  data  points  are  not 
obscured.  For  the  time-of-flight  measurements,  the  error  bars  represent 
sound  speed  uncertainty  due  to  uncertainty  in  the  sensor  separation  distance 
and  the  time  resolution  of  the  data  acquisition.  The  EDFM  was  evaluated 
using  the  parameters  in  Table  II  and  the  mean  values  of  porosity  and  per¬ 
meability. 


sound  speed  are  repeated  here  for  convenience.  The  EDPM 
considers  a  water-saturated  granular  sediment  with  a  total 
mass  density  p,  pore  fluid  mass  density  p;,  sediment  particle 
mass  density  p5,  individual  sediment  grain  bulk  modulus  Kn 
pore  fluid  bulk  modulus  Kp  porosity  /3,  tortuosity  ay  perme¬ 
ability  ac,  pore  fluid  viscosity  77,  and  complex  sound  speed  c. 
The  sound  speed  is  a  function  of  angular  frequency  co  and  is 
given  by 


c  = 


H 


PeffM  ’ 


(3) 


where  the  compressibility  is  determined  by  a  mixture  rule 


H  = 


1  zl 

Kr 


K, 


and  the  effective  density  is 


(4) 


TABLE  II.  Material  parameters  used  in  the  effective  density  fluid  model  (SI 
units).  Parameters  obtain  from  measurements  on  the  sediment  are  in 
BOLDFACE.  The  mean  value  is  followed  by  the  range  of  variation  among 
multiple  measurements.  The  remaining  parameters  are  from  Ref.  12. 


Porosity 

0.370  (0.368-0.371) 

Sand  grain  mass  density 

2650 

Water  mass  density 

998 

Sand  grain  bulk  modulus 

3.6  X1010 

Water  sound  speed 

1492 

Water  viscosity 

0.001 

Permeability 

4.58X  10"11  (4.18  X  10_11-4.74X  10"n) 

Tortuosity 

1.25 
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TABLE  III.  Material  parameters  used  in  the  evaluation  of  Eq.  (Al)  (SI 
units). 


Sediment  mass  density,  pL 

2039 

HDPE  mass  density,  pw 

1000 

HDPE  compressive  sound  speed,  c{ 

2500 

HDPE  transverse  sound  speed,  c, 

680 

/ 


PefM  =  p, 


a{\  -  P)ps  +  P(a-  1  )pf+ 


ifipFy 

pf<i)K 


1  -  p)ps  +  (a  -  2/3  +  + 


ipFrj 


(5) 


The  parameter  F  represents  a  deviation  from  Poiseulle  flow 
of  the  fluid  within  the  pores  and  is  given  by 


F(e)  = 


;  He) 

4 

2T“- 

1  -  -T(e) 
€ 


(6) 


where  the  Bessel  functions  Jv(x)  are  used  in 

,  (-  \<Vi(eV0 

7(e)  = - ^ - , 

Jo(c'  0 


(7) 


where  the  quantity  e=a\o)pj/ 77  in  conjunction  with  <2 
=  \8 a/c/yS  represents  the  pore  size  parameter.  The  sound 
speed  measurements  shown  in  Fig.  6  are  compared  to  the 
model  phase  speed  given  by  real  part  of  Eq.  (3)  and  the 
material  parameters  of  Table  II. 


VI.  EFFECT  OF  FINITE  WALL  IMPEDANCE 

An  exact  analytic  model21  for  sound  propagation  in  an 
elastic-walled,  fluid-filled  cylindrical  tube  was  used  to  inves¬ 
tigate  the  effects  of  finite  wall  impedance  on  the  sound 
speeds  inferred  from  the  resonance  frequency  measurements 
for  the  axisymmetric  modes.  The  model  yields  particle  dis¬ 
placement  field  equations  and  a  dispersion  relation  for  the 
axisymmetric  modes  of  an  inviscid-liquid-filled  cylindrical 
tube  with  arbitrary -thickness  elastic  walls,  where  the  outer¬ 
most  radial  boundary  condition  was  approximated  as  pres¬ 
sure  release.  The  dispersion  relation  has  been  verified  for  a 
variety  of  sizes  and  wall  materials,  including  steel,17  alumi¬ 
num,  and  PVC,22  and  is  presented  in  Eq.  (Al). 

The  unconfined  (intrinsic)  sound  speed  of  the  material 
that  fills  the  tube  is  an  input  parameter  to  Eq.  (Al)  and  the 
admissible  phase  speeds  are  calculated.  This  model  is  limited 
to  axisymmetric  modes,  so  we  applied  it  to  the  0, 1  ,N  modes 
and  found  that  the  maximum  deviation  of  the  sediment  sound 
speed  inferred  from  the  pressure-release  model  was  within 
0.57%  of  the  sound  speed  inferred  from  the  elastic  wave¬ 
guide  model  and  the  minimum  deviation  was  0.18%.  The 
material  parameters  used  for  the  evaluation  of  Eq.  (Al)  are 
given  in  Table  III  and  the  resulting  sediment  sound  speeds 
are  given  in  Table  IV.  The  corresponding  sediment  sound 
speeds  obtained  with  the  pressure-release  model  are  also  in¬ 
cluded  in  Table  IV  for  comparison. 


TABLE  IV.  Comparison  between  sediment  sound  speed  values  c  of  Fig.  6, 
which  were  inferred  from  the  pressure -re lease  model  [Eq.  (2)],  and  cor¬ 
rected  sound  speeds  c]t  obtained  from  the  axisymmetric  elastic  waveguide 
model  [Eq.  (Al)]. 


Frequency  (Hz) 

c  (m/s),  Eq.  (2) 

c]  (m/s),  Eq.  (Al) 

2179 

Circle  (°)  data  points 
1721 

1712 

2892 

1736 

1731 

3788 

1742 

1739 

2178 

Square  (□)  data  points 
1721 

1711 

2899 

1741 

1735 

3781 

1739 

1736 

2186 

Cross  (X)  data  points 
1727 

1717 

2899 

1741 

1735 

2160 

Diamond  ( O )  data  points 
1706 

1697 

2865 

1720 

1714 

The  relative  size  of  the  measurement  uncertainty  due  to 
tank  dimension  uncertainty  is  compared  to  the  finite  wall 
impedance  correction  in  Fig.  7.  For  clarity,  this  comparison 
is  only  shown  for  one  of  the  data  sets.  The  circle  data  points 
from  Fig.  6  are  repeated  in  Fig.  7  and  their  values  after 
correction  for  the  finite  wall  impedance  are  shown.  The  finite 
wall  impedance  biases  the  axisymmetric  resonance  measure¬ 
ments  (circles  and  squares  in  Fig.  6)  by  +0.57%  or  less,  but 
this  is  significantly  less  than  the  measurement  uncertainty 
due  to  the  tank  dimensions. 

An  exact  analytic  model  for  asymmetric  excitation  of  a 
liquid-filled  waveguide  with  walls  of  arbitrary  thickness  is 
not  available.  We  used  an  approximate  model^  which  treats 
the  wall  as  a  shell  and  yields  a  dispersion  relation  Eq.  (A6) 
for  both  axisymmetric  and  asymmetric  modes.  Otherwise, 
the  corrections  procedure  is  the  same  as  for  the  axisymmetric 
case,  discussed  earlier.  The  dimensions  of  the  sediment  tank 


FIG.  7.  The  effect  of  ignoring  the  finite  impedance  of  the  tank  walls  is 
shown  for  the  axisymmetric  modes.  The  large  circles  are  the  sediment  sound 
speeds  inferred  from  a  pressure -release  model  of  the  system,  repeated  here 
from  Fig.  6,  and  the  error  bars  represent  the  measurement  uncertainty  due  to 
imperfect  knowledge  of  the  tank  dimensions.  The  small  circles  are  the  sedi¬ 
ment  sound  speeds  inferred  from  a  model  which  accounts  for  the  finite 
impedance  of  the  walls  of  the  tank.  Ignoring  the  finite  impedance  of  the 
walls  results  in  a  maximum  bias  error  of  +0.53%  at  the  lowest  frequency 
measured. 
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FIG.  8.  The  data  and  model  prediction  shown  in  Fig.  6  are  repeated  here, 
except  that  the  resonance  method  measurements  were  corrected  for  the  finite 
wall  impedance  of  the  tank.  The  extent  of  measurement  uncertainty  for  the 
resonance-based  measurements  (due  to  sensor  position  and  tank  geometry 
uncertainty)  is  indicated  by  the  closed  polygon  surrounding  the  data  points. 
The  solid  curve  is  the  EDFM  prediction  based  on  the  mean  porosity  and 
mean  permeability.  The  dashed  curves  above  and  below  the  solid  curve 
represent  the  maximum  and  minimum  EDFM  predictions  associated  with 
the  range  of  the  porosity  measurements  and  the  standard  deviation  of  the 
permeability  measurements. 

are  well  within  the  limits  required  by  the  shell  approxima¬ 
tion.  To  ensure  continuity,  we  compared  the  phase  speeds 
calculated  by  Eqs.  (Al)  and  (A6)  for  the  symmetric  mode 
0,1  and  found  excellent  agreement.  We  are  therefore  confi¬ 
dent  that  the  approximate  Eq.  (A6)  is  accurate  for  the  geo¬ 
metric  and  material  properties  used  here.  The  data  and  model 
prediction  presented  in  Fig.  6  are  repeated  in  Fig.  8,  but  with 
the  finite  wall  impedance  correction  applied  to  all  the  data 
from  both  axisymmetric  and  asymmetric  modes.  The  asym¬ 
metric  data  before  and  after  correction  are  presented  in  Table 
V,  and  shows  a  maximum  correction  of  +0.76%.  Note  that 
the  time-of-flight  measurements  are  not  biased  by  the  pres¬ 
ence  of  the  tank  walls  because  the  environment  appears  to  be 
a  free  field  for  the  pulse  lengths  and  frequencies  used. 

VII.  CONCLUSIONS 

Laboratory  sound  speed  measurements  were  obtained  in 
a  water-saturated  granular  sediment  composed  of  commer¬ 
cial  silica  blasting  sand  and  fresh  water.  Above  20  kHz,  a 
direct  time-of-flight  method  was  used  and  between  2  and 
5  kHz,  an  indirect  method  was  used  to  infer  frequency- 
dependent  sound  speeds  from  measured  resonance  frequen- 

TABLE  V.  Comparison  between  sediment  sound  speed  values  c  of  Fig.  6, 
which  were  inferred  from  the  pressure-release  model  [Eq.  (2)],  and  cor¬ 
rected  sound  speeds  Cj,  obtained  from  the  axisymmetric  elastic  waveguide 
model  [Eq.  (A6)]. 


Frequency  (Hz) 

c  (m/s),  Eq,  (2) 

C\  (m/s),  Eq.  (A6) 

3207 

Cross  ( X )  data  points 
1723 

1710 

3747 

1740 

1730 

3195 

Diamond  ( 0 )  data  points 
1716 

1704 

3735 

1734 

1725 

cies.  The  pressure-release  approximation  yields  a  surpris¬ 
ingly  good  model  for  a  thin-walled  plastic  tank  filled  with  a 
water  saturated  sediment,  as  shown  in  Figs.  6-8.  There  is  at 
maximum  about  0.5%  error  in  sound  speeds  inferred  from 
the  pressure  release  model. 

Considering  the  more  realistic  analysis  which  does  ac¬ 
count  for  the  finite  impedance  of  the  tank  walls  and  is  shown 
in  Fig.  8,  the  EDFM  model  “  describes  the  sound  speed  mea¬ 
surements  very  well.  Most  of  the  low-frequency  measure¬ 
ments  obtained  with  the  resonance  method  are  just  lower 
than  the  predicted  values,  but  the  agreement  between  model 
and  measurement  is  well  within  the  measurement  uncer¬ 
tainty.  Three  of  the  four  high-frequency  time-of-flight  mea¬ 
surements  are  just  higher  than  the  prediction,  but  again, 
within  the  range  of  measurement  uncertainty.  The  overall 
conclusion  is  that  dispersion  is  clearly  observed,  as  predicted 
by  the  Biot-derived  models,  and  that  the  EDFM  accurately 
describes  the  sound  speed  observed  in  this  artificial  labora¬ 
tory  sediment.  Note  that  within  the  experimental  frequency 
range,  both  measurements  and  model  indicate  about  4%  dis¬ 
persion.  Evaluation  of  the  EDFM  from  its  low  frequency 
asymptote  (—100  Hz)  to  its  high  frequency  asymptote 
(— 1  MHz)  yields  a  total  of  about  8%  dispersion.  A  second¬ 
ary  conclusion  can  be  drawn  from  the  spread  in  the 
resonance-inferred  data,  which  sampled  the  entire  sediment 
volume,  as  compared  to  the  high  frequency  measurements 
which  only  sampled  a  small  portion  of  the  volume.  Even  in  a 
carefully  prepared  laboratory  specimen,  granular  sediments 
are  acoustically  nonhomogeneous. 

Overprediction  of  sediment  sound  speed  by  Biot-type 
models  in  the  transition  region  between  the  low  and  high 
frequency  asymptotes  has  been  previously  reported  in  Refs. 
1  and  4.  One  explanation  for  this  deviation  in  the  previous 
work  is  the  existence  of  undissolved  air.  In  the  present  ex¬ 
periment,  significant  effort  was  expended  to  achieve  a  gas- 
free  sediment,  and  it  appears  that  if  there  was  any  remaining 
undissolved  air  (in  microscopic  crevices  for  instance)  the  ef¬ 
fect  was  not  acoustically  significant. 
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APPENDIX 

1.  Axisymmetric  dispersion  relation 

We  considered  sound  propagation  along  the  longitudinal 
axis  of  a  lossless  fluid  cylinder  enclosed  by  a  lossless,  arbi¬ 
trary  thickness,  elastic  cylindrical  wall.  Acoustic  energy  at 
angular  frequency  o)= 2 irf  propagating  in  this  system  (where 
/  is  the  number  of  cycles  per  second)  is  governed  by  a  set  of 
exact  wave  theoretic  equations  given  originally  by  Del 
Grosso"1  and  refined  by  Lafleur  and  Shields.22  The  disper¬ 
sion  relation  for  this  system  is  given  in  Eq.  (Al).  The  ele- 
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ments  of  Eq.  (Al)  can  be  complex  valued,  but  real  values  of 
qom  =  (ol cQm  that  cause  the  left-hand  side  of  Eq.  (Al)  to  be¬ 
come  zero  lead  to  propagating  modes  that  have  phase  speed 
c0m.  At  any  given  frequency  there  can  be  multiple  permis¬ 
sible  modes,  which  are  indexed  by  integer  values  of  m.  Since 


the  phase  speed  is  frequency-dependent,  at  a  given  sediment 
cylinder  resonance  frequency  /0>wA/,  the  phase  speed  will  at¬ 
tain  a  value  denoted  by  c0ymN.  The  zero  subscript  that  pre¬ 
cedes  the  index  m  indicates  that  these  are  axisymmetric 
modes.  The  dispersion  relation  is 


1  +  [Lu(Pm)L00(Tm)\ 
/  irbdPlTi 


8  El 


+  [Z.„(7m)L00(Pm)] 


l^rbdE2m 


+  [L]0(Pm)L0](Tm)  +  L0l(Pm)Lw(Tm)] 


8 


+  d(\+Qmb)Ln(Tm)L0l(Pm)] 


+  [bLu{Pm)Lw(Tm)  +  d(  1  +  Qmb)Lu(Pm)L0i(TJ] 
rrPmEm  7 rP 


t rP2mTm  t rPiqlmT„ 

m  rri  _ m~\)m 

8  Em 


8  E~ 
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8 


+  [(l+emWLn(7^)L11(|>m)] 


irPm  irpjql  _  7 
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=  0, 
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r 


where  the  following  definitions  have  been  used: 


*om 

■  Pm 

=  v*?  -  q\m.  Tm  = 

\k2- 

<7cL- 

(A2) 

Q0m 

=  (o/c0m,  kx 

=  w/c'i 

,  ki  -  col  ch  kt  - 

c ole ,, 

(A3) 

Em  ' 

=  q20m-k2/2 , 

Qm  = 

pLu>2bJ0(X0m) 

(A4) 

2pwC'X0mJt(X0m) 

Emn 

(0  =  Jm(d$Y, 

■Jn(b&Ym(d&, 

(A5) 

and  the  geometric  and  material  input  parameters  are:  c{  is  the 
intrinsic  sound  speed  in  the  inner  liquid  (the  sound  speed  that 
the  sediment  would  exhibit  in  the  free  field),  c/  and  ct  are  the 
longitudinal  and  transverse  (shear)  speed  for  the  wall  mate¬ 
rial,  pi  and  pw  are  the  liquid  and  wall  material  densities,  and 
b  and  d  are  the  inner  and  outer  radii  of  the  pipe,  respectively. 
yp(jc)  and  Yu(x)  are  Bessell  functions  of  the  first  and  second 
kind,  respectively. 

2.  Asymmetric  dispersion  relation 

We  modeled  nonaxisymmetric  sound  propagation  along 
the  longitudinal  axis  of  a  lossless  fluid  cylinder  enclosed  by 
a  lossless,  finite  thickness,  shell-like  cylindrical  wall,  using 
the  formulation  of  Fuller  and  Fahy.2^  The  dispersion  relation 
is  given  in  terms  of  a  determinant 

det[L,j]  =  0  (A6) 

where 

Ln  =  -n2  +  (A'n»2  +  5(l  -  v)n 2, 

Ln=\{\  +  v)n(knmb),  Ln=v(knmb),  (A7a) 

L2i  =  L12,  L12=-Cl2 +  \{\  -  v){knmb)2  + n2,  L22  =  n, 

(A7b) 


—  T13,  />32  —  £03  » 

L33  =  -  il2  +  1  +  f?[(knmb)2  +  n2]2  -  FL.  (A7c) 

In  Eq.  (A7),  11  is  the  dimensionless  frequency,  Cl  =  cob/ch  v 
is  the  ratio  of  Poisson  of  the  wall  material,  /3  is  the  wall 
thickness  parameter  02=h2/  \  2b2,  h  is  the  wall  thickness,  and 
FL  represents  fluid  loading 

FL  =  Cl2(pL/pw)(b/h)(  1  /krmb)[Jn(krmb)/J'n(krmb)] .  ( A8) 

The  radial  fluid  wave  number  ksmb  and  the  axial  wave  num¬ 
ber  knmb  are  related  by 

krmb  =  [(l2(c,/cl)2-(knmb)2]m.  (A9) 


3.  Illustration  of  finite  wall  impedance  correction 
scheme 

The  effect  of  the  finite  wall  impedance  and  the  method 
of  correction  for  this  effect  are  demonstrated  in  Fig.  9  for  the 
axisymmetric  mode.  The  phase  speed  for  the  0, 1  ,jV  mode  is 
shown  for  both  the  pressure  release  model  and  the  elastic 
waveguide  model  Eq.  (Al).  Note  that  the  predicted  phase 
speeds  are  similar,  therefore  we  do  not  expect  the  corrections 
to  be  large.  For  a  given  experimentally  observed  resonance 
frequency  /0  1>Ar,  the  observed  modal  phase  speed  is  calcu¬ 
lated  from 

^o,\,n=  \Zo.ijv»  (A  10) 

where  \=2 h/N.  Then  the  effective  sediment  sound  speed  cx 
is  varied  in  Eq.  (Al)  until  the  predicted  modal  phase  speed  is 
Cq.\,n-  This  procedure  was  repeated  for  each  resonance  fre¬ 
quency  identified  from  an  axisymmetric  mode.  The  correc¬ 
tion  procedure  is  identical  for  the  asymmetric  modes,  except 
Eq.  (Al)  is  replaced  by  Eq.  (A6).  The  results  are  reported  in 
Table  IV  and  shown  in  Fig.  8. 
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(A  12) 


X  =  2 h  X  =  2/i/3 

/0iU  =2179  Hz  /0ili2  =  2892  Hz  /o.u  =  3788  Hz 


FIG.  9.  The  effect  of  the  finite  wall  impedance  on  the  modal  phase  speed 
Cq.ijv-  The  vertical  axis  is  normalized  by  the  sediment  effective  intrinsic 
sound  speed  Cj.  The  basis  for  the  procedure  used  to  correct  for  the  finite  wall 
impedance  is  illustrated  for  the  one  of  the  data  sets  (the  open  circles  of  Fig. 
6).  Integer  multiples  of  one-half  wavelength  occupy  the  sediment  cylinder 
along  its  height  for  the  0,1, N  modes.  Knowledge  of  these  mode  shapes 
from  Eq.  (1),  the  measured  resonance  frequencies  /o.i^v  and  Eq.  (A  10)  allow 
one  to  use  Eq.  (Al)  to  relate  the  measured  value  of  c0)ijv  t0  intrinsic 
effective  sediment  sound  speed  c 1  that  would  be  measured  in  the  free  field. 
The  results  of  the  correction  are  shown  in  Table  IV 


4.  Effect  of  neglecting  viscosity  in  the  wall 
impedance  correction 

The  above-presented  elastic  waveguide  models  approxi¬ 
mate  the  internal  fluid  as  inviscid,  therefore  the  fluid  can  slip 
at  the  boundary  with  the  wall.  In  this  section  it  is  shown  that 
neglecting  the  sediment  viscosity  in  the  waveguide  models 
has  an  insignificant  effect  on  the  phase  speeds  in  this  expen- 
ment.  This  is  due  to  the  large  size  of  the  sediment  tank  in 
comparison  to  the  boundary  layer  at  the  wall. 

In  this  work,  the  sediment  is  approximated  as  an  effec¬ 
tive  fluid  and  can  be  assigned  an  effective  viscosity  with 
respect  to  pipe  wall  interaction,  as  is  done  in  the  modeling  of 
slurry  flows,  4  and  in  acoustic  studies  of  sediments  in 
tubes. 25  For  mixtures  of  liquids  and  sand  with  high  solid 
volume  concentrations  (such  as  the  present  sediment)  the  ef¬ 
fective  mixture  viscosity  jUmix  is  dependent  upon  the  solid 
volume  concentration  <£  as 


—  =  1+2.50  +10.05^  +  2.73  X  l(T3exp(  16.60),  (All) 

where  fi  is  the  viscosity  of  the  host  liquid/4  The  sediment 
porosity  in  this  work  is  0.37  therefore  <£=0.63  and  Eq. 
(All)  yields  a  viscosity  ratio  Mmix//x=  102. 

For  acoustic  waves  propagating  within  a  viscous  fluid 
inside  a  pipe,  an  acoustic  boundary  layer  is  formed  at  the 
wall.  The  resulting  acoustic  effects  can  be  modeled  with  a 
complex  wave  number"6 


i )  ^wall » 


where  c  is  the  effective  sound  speed  of  the  sediment  inside 
the  waveguide  and  cu/c>awal \  is  assumed.  Thermal  losses 
are  considered  negligible  compared  to  the  viscous  losses, 
which  are  characterized  by 

_  J_  Mmix^  /  a  i 

"wall-fl\  2pc2  ’  3 


where  a  is  the  inner  radius  of  the  pipe  and  p  is  the  total  mass 
density  of  the  sediment.  Equation  (A  12)  indicates  a  reduced 
phase  speed  inside  the  waveguide  due  to  the  viscous  bound¬ 
ary  layer.  The  effect  of  the  viscosity  on  phase  speed  can  be 
expressed  in  a  ratio 


£  _  /  i  _j_  ££Wn\ 

c  \  0)  / 


(A  14) 


where  c-Rc[a)/k]  is  normalized  by  the  sound  speed  in 
absence  of  a  viscous  boundary  layer,  c. 

Using  the  sediment  mass  density  from  Table  III,  the 
highest  frequency  in  the  resonance  regime  of  the  experiment 
(3788  Hz)  and  the  mean  sound  speed  observed  inside  the 
waveguide  at  the  same  frequency  (1740  m/s),  Eq.  (A  14) 
yields  c/c- 0.99905,  which  is  a  negligible  systematic  error  in 
comparison  to  other  sources  of  measurement  uncertainty  in 
this  experiment  and  is  ignored. 
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